Spin-wave softening and Hund's coupling in ferromagnetic manganites 
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Using one-orbital model of hole-doped manganites, we show with the help of Holstein-Primakov 
transformation that finite Hund's coupling is responsible for the spin-wave softening in the ferromag- 
netic -B-phase manganites. We obtain an analytical result for the spin-wave spectrum for Jh 3> t. 
In the limit of infinte Hund's coupling, the spectrum is the conventional nearest-neighbor Heisen- 
berg ferromagnetic spin- wave. The o(t / Jn)-order correction is negative and thus accounts for the 
softening near the zone boundary. 

PACS numbers: 75.30.Vn,71.10.-w 
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The observations of large magnetoresistence (LMR) 
in Ndo.5Pbo.5Mn03 , giant magnetoresistence (GMR) 
and colossal magnetoresistence(CMR) in manganites 
(Ri_ a A x Mn03, R is a rare earth element and A a di- 
valent alkaline-earth metal) a decade ago pjhave rekin- 
dled much interest in these materials which have been 
known for half a century [2j. Upon doping, the mangan- 
ites undergo complicated transitions resulting in various 
magnetic, charge-ordering and orbital-ordering phases, 
showing the interplay between relevant spin, charge and 
orbital degrees of freedom. In particular, magnetism 
and electronic transport are clearly correlated. So it 
is widely believed that knowledge of spin dynamics can 
provide important information of the underlying physics 
of CMR. Perring et al first measured the spin waves in 
Lao.7Pbo.3MnC>3 for a broad range of g[3]. The magnon 
spectrum is well defined at low temperatures and can be 
accounted for by the nearest neighbor Heisenberg model. 
Subsequent measurements for Pr .63Sr .37MnO3 and 
Nd o .7Sr o . 3 MnO30], Ndo.7Ba .3Mn0 3 Q further showed 
that the magnon spectrum deviates from the Heisenberg 
model and becomes softened near the zone boundary. So 
the behavior seems a universal phenomenon of mangan- 
ites. 

As is well known, a number of interactions such as 
spin-orbital coupling, Hund's coupling, antiferromagnetic 
coupling between core spins, Coulomb interaction and 
dynamic Jahn- Teller effect coexist in manganites. These 
interactions are supposed to explain the existence of dif- 
ferent phases of doped manganites. To explain the spin 
wave softening, various mechanisms were proposed. The 



authors of Q further showed that the experimental spec- 
trum can be reproduced reasonably well by an extended 
Heisenberg model. Furukawa argued that the soft- 
ening seems to be explainable by ferromagnetic Kondo 
lattice model with bandwidth narrower than the Hund's 
coupling. Solovyev et al 7] showed that the spin-wave 
behavior near the zone boundary has a purely magnetic 
spin origin, and neither the lattice deformation nor the 
orbital ordering are required to account for the soften- 
ing. Dai et al argued that the observed magnon soft- 
ening and broadening are due to strong magnetoelastic 
interactions 8J . And this magnon-phonon coupling was 
later treated quantitatively inj^]. Using ferromagnetic 
Kondo lattice model and composite operator method, 
Mancini et al obtained the softening spectrum 10] . Shan- 
non et al constructed a theory of spin wave excitations 
in the bilayer manganite Lai.2Sri.sMn207 based on the 
simplest double exchange model and explained partly the 
softening behavior Krivenko et al showed that the 
scattering of spin excitations by low-lying orbital modes 
may cause the magnon softening 12]. 

In this paper, we show that in the hole-doped man- 
ganites, the softening behavior might be of a purely elec- 
tronic origin ,i.e., a strong but finite Hund's coupling 
between the e g electron and the core spin. Since in the 
hole-doped manganite there is less than one e g electron 
per site on average and the d x 2_ y 2 orbital energy is signif- 
icantly higher than that of d 3z 2_ r 2 |l3| due to Jahn- Teller 
splitting, one orbital description is a reasonable approx- 
imation. As in |14j . we adopt the model Hamiltonian 



a i <i,j> iff i 
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where t is the double exchange hopping, are near- represents the e g electrons, Jh is the Hund's coupling 
est sites, ju is the chemical potential for the fermions, c icr between the e g spin s ; = \c\crci and the the core spin 
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Si- Jaf is the antiferromagnetic interaction between the 
core spins, which is necessary to account for the G-phase 
parent (x = 1) manganites. The last term is the Hubbard 
Coulomb interaction. We use the Holstein-Primakoff 
transformation for the core spins (S = 3/2); 5j + = 
(2S-ala i ) 1 / 2 a i ,Sr = ajf^S-aJoi) 1 /*, S? = S-atai and 



take the approximation {2S-a\a i ) 1 l 2 ~ (2S 1 - (ajai)) 1 / 2 . 
Homogeneity implies that {a\a,i) = (a^a). Because we 
consider the low temperature case, we can drop the 
magnon quadratic term, hence the total Hamiltonian can 
be written as 



H = <^^cLcj ff -/i^4 £7 c iff + C/^n iT n u -ij H v4^(s i f a 1 t + Si ai )-J H S^s? 

a \a i i i 

+A 2 J AF a;a j + JafZNS 2 - 2ZJ AF S ajas + J H ^ s-ajai 

(U> i i 



(2) 



, where ^4 2 = 2S* — (a^a), Z = 6 is the coordination num- method, we consider the doublet B(i) = (cti, sj l ~) T . The 
ber of the core spins. To use the composite operator equation of motion for B(i) is 



id t B(i) = [B(i),H] 



-i JnAsf + A 2 J AF E e a i+e - 2ZJ AF S ai + J H s? ai 



t Ee( C iT C i+el - c i+eT C u) ~ J H As i a i + JkSs{ - J H sta[ai 
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(3) 



Composite operator method assumes that the right-hand 
side can be expressed as 



[S(i),#] = £>(i,j)B(j) 
j 

with e(i, j) determined in the following way, 
£ (i,j)=^m(i,l)/- 1 (l,j) 



(4) 
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where J(i,j) = <[B(i),flt(j)]) )TO (ij) 
([ic?t-B(i), and (} represents the expectation 

value. Thus e(i,j) contains some parameters to be de- 
termined self-consistently. This approach was proposed 
for Hubbard model originally and recent intensive 
studies0| show credible agreement with Monte Carlo 
method. In our case (again due to homogeneity, 

J(i,j) = <5 ij -diag(l,2( s z )) 
m n (i,j) = ^(Jh^ 2 ) -2ZSJ AF ) + A 2 J AF Y,Sj, i+e 

e 

mi 2 (i,j) = <5ijJ H (-A(s 2 ) - (sr ai )) 

m 22 (i, j) = -tpi - <?>j,i+e) + JnA(s7ai) 

e 

+2J H 5( S 2 ) -2J H ( S ?a^i) 



where pi = J2a( c L^) ^ P^ = ( s i a i)>Ps = (sfaja ; }. In 
the k-space 

m n (k) = (J H (s 2 )-2Z,SJAF) + ^ 2 JAF7k 
mi 2 (k) = J H (-A(s 2 ) -p 2 ) 

m 22 (k) = J n Ap 2 -tZp 1 (l-j hL ) + 2J Fl S(s z } -2JhP 3 

We assume that at T = OK, (a^a) = 0, which satisfies self- 
consistency using the resulting retarded Green's function 
and spectral theorem. Then condition W|k=o = requires 
that P3 = — \p2{A + 7§§y)' So the e-matrix is 

eu (k) = J H (s 2 ) -2Z5J A F(l-7k) 

£l2(k) = -^(A{ S *)+p 2 ) 

eai (k) = -J H (A( s 2 )+p 2 ) 

£ 22 (k) = J H 5-^(l- 7 k) + JH^7^V + J H 



and the Green's function is 

w - (JuS - |^(1 - 7k) + JhA^ + J H 



Dn(w,k) = 
Di 2 (w,k) = D 2 i(k) = 
£ 22 (w,k) 



p 2 



2 <s 2 => 2 



(w - cJi(k))(w - w 2 (k)) 
-J H (A(s 2 )+ P2 ) 
(a; - wi(k))(o; - w 2 (k)) 
2(s z )[u; - (J H (s 2 ) - 2Z5J AF (1 - 7k)) 



(w - wi(k))(w - w 2 (k)) 
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where wi^k) are acoustical and optical branches of the 
spin excitations. Using 



1 x - t 
A? 2-^, 
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£>i2 (w + if], k) - D12 (w - irj, k) 



(6) 

we have at T = OK, p 2 = 0, therefore, p 3 — 0. Ac- 
cordingly, in this scheme, there are two parameters left: 
(s z ),pi. The acoustical magnon spectrum can be ex- 
panded as a Taylor series which manifests the role of 
Hund's coupling 



Wl (k)=t]T(^)"a„(l- 7k ) 



n+l 



n=0 



Jh' 



(7) 



where as usual, 7k = Z 1 ^ £ 
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The first few a n are 
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(s z ) + 5 
9S(4S^(s z ) - pi) 2 
<s z )((s*) + 5) 3 
27(4S^(s z ) - Pl ) s S(S - (s z )) 



a z\2 



8l(S 2 - 3S{s z ) + {s z ) 2 )(4S^(s z ) - Pl fS 



(s z ) 3 (S + (s z )) 7 



In the small k limit, 



uj ~ Dk 2 ,D = -(45 2 ^ + 
pi)/((s z ) + 5). Note that the hopping energy tp\ is neg- 
ative and when it overcomes the AF term, the resulting 
magnon stiffness D is positive. Our numerical results 
show that this self-consistency is satisfied. Expression 
(7) suggests that the softening comes from the finite Jh- 
To fix the parameters (s 2 ), we use spectral theorem and 
get (s z ) — 4(1 — x), where x is the dopant concentration. 
To fixpi, we need the fermion sector. Using the notations 
in [ljj f° r t ne fermion operator ip(i) = 77-pi, ?7xi) T i 
where £ CT = (1 — n._ ff )c CT , r\ a = n^ a c a are the Hubbard 
operators, we obtain the retarded Green's function for ip 
in the large— U limit at zero temperature. 



G R (u,k) =diag( 



,0, 



(8) 

('diag' means diagonal matrix) with -Ei(k) = — /i+6i 7k — 
|5J Hj -E 2 (k) =U - fi + 6tu + 6to 7k , E 3 (k) = [24ta 7k - 
2fix + SJ H x + 12t Pl j k - 12t 7k + 12iA T ]/(2x), E 4 (k) = 
[-2Ux + 2[ix - SJ R x + 2U + 12ip i7k + 5J H - 2fi + 
12fAj]/[2(l — x)] , where A is related to the nearest- 
neighbor correlations of the Hubbard operators : A a = 
(^ <T (i+e)^+(i))-(i7 (r (i+e)r/t(i)). In this scenario, E 1} E 3 
are partially filled and E2 , E4 are empty. The relevant 
parameters are [i, A-f,pj_. We have three equations to 
fix them 1 - x = 2 - Cf x - C& - Cf 3 - Cf 4 , , A T = 
= C 3 F 3 , where C F = WW(i)), = 
{ip(i + e)^t(i)>. We know that C[ 2 = and Cf 4 = 1 - x. 



Thus Cf x = x = C33, so -E3 is empty, i.e., only £1 is 
partially filled. Hence only /i is relevant to our problem 
and can be fixed by x — N^ 1 ^2 k 0(Ei(k)) where 9{x) is 
the usual step function. The hopping energy is 

t x - 



tpi 



F 7 



N ^ 

k 



- J B 1 (k)) 7k <Q 



(9) 



The other two parameters A|,p| can also be determined 
by iA T = -tpi > 0, 24tx + 12t( Pl - 1) = -2^x + 5 J H x + 
12tAf. Further analysis show that for Jh > 2.5t, the 
whole scheme is self-consistent. Fig.l shows that two rel- 
evant fermion bands for x — 0.301, t — 1, Jh = 3.0 ( in 
unit of t). 

It is seen from the magnon spectrum (7) that we can 
estimate the two model parameters t and Jh from mea- 
sured data. Fig. 2 shows the comparison between our cal- 
culated result for the prescribed antiferromagnetic cou- 
pling Jaf = 0.01 and the measured result at T = 10K 
for Pr .63Sr .37MnO3 in 4). The solid curve in the left 
panel is the fit to a nearest -neighbor Heisenberg model 
and gives the value at zone boundary about 34.2meV. 
This corresponds to the the uppermost curve in the 
right panel. The comparison gives the hopping energy 
t ~ 0.462eV. The circles are the data measured and give 
the value at zone boundary about 23meV, corresponding 
to the point 0.05 in the right panel. This point corre- 
sponds to Jh — = 3.2t ~ 1.48 eV. Note that the ratio 
Jn/t is very close to the values of interaction from a 
number of references (lr| (l7| . It is worth noting that the 
nearest-neighbor Heisenberg interaction alone can not ac- 
count for the Curie temperature. The fitting curve in the 
left panel corresponds to the nearest-neighbor Heisenberg 
spectrum w(k) ~ 51.3(1 — 7k ) meV. In the mean field the- 
ory,the Curie temperature T c corresponding to the spec- 
trum w k = 2ZS'*J(1 - 7k ) is k B T c = pZS*{S* + 1) 
(here S* = S+ |(1 — x) is the effective spin). This gives 
T^ IF ~ 500K. Taking into account that in three dimen- 
sions for a simple cubic lattice , the real Curie temper- 
ature T c and T^ F have relation^: T c = 0.75T«¥ F , 
we get T c ~ 375 K, higher than the real value 315 K. 

To conclude this paper, we present some discussions 
and comments. In the derivation of the series expres- 
sion of the magnon spectrum, we have used the ap- 
proximation (25 - alai) 1 / 2 ~ (25 - (atai)) 1 / 2 in the 
Holstein-Primarkov transformation. This can be sat- 
isfied at very low temperatures. Further, the quartic 
term Jaf E(i j) °i a i a j a j is neglected because Jaf is very 
small and the magnon fluctuation at zero temperature 
is negligible. The series expression (7) of the accous- 
tic magon dispersion shows alternating behavior; conver- 
gence is guaranteed when Jh5 > 3. The model param- 
eters t, Jh, Jaf and the hopping energy pi can be esti- 
mated by fitting experimental data. There is a simple 
physical picture for the deviation of the magnon spec- 
trum from that of Heisenberg model. The interaction be- 
tween core spins is induced by the hopping of e g -electrons 
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FIG. 1: The bands for the Hubbard operators at T = OK. 
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behavior is still an issue of debate. Based on the ob- 
served proximity of phonon dispersion and magnon dis- 
persion and the anisotropic spin-wave broadening, Dai 
et alQ argue that strong magnon-phonon coupling is 
needed for a complete understanding of the low tempera- 
ture spin dynamics of manganites. Quite recently, Endoh 
et al concluded 20] that the ferromagnetic magnons in 
Smo.55Sro.45Mn03 is of orbital nature since the magnon 
dispersion shows anisotropy which is mainly determined 
by the short range correlation of the e g orbitals. They 
explained that the anisotropic magnon dispersion are at- 
tributed to long range magnetic interactions based on 
fitting the data to Heisenberg model with long range 
interactions. We believe that if orbital degrees of free- 
dom are taken into account in our model, the resulting 
magnon spectrum will be anisotropic since orbital degrees 
of freedom bring anisotropy into the system. Finally, we 
remark that as manganites are very complex systems, 
there might be multiple mechanisms contributing to a 
single phenomenon. The analysis provided in this paper 
shows that Hund coupling might be of primary impor- 
tance. 
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FIG. 2: The spin wave spectrum along FX of _B-phase. The 
left panel (taken fromQ) is the experimental result. The right 
panel is the calculated spin-wave spectrum (in unit of t) at 
T = OK . 

and the dominant term is linear in t. If the Hund's cou- 
pling Jh is infinite, only the dominant term plays the role. 
The e g electron and core spin must add up to a total spin- 
2 to minimize the energy in the -B-phase. So the actual 
background for spin excitation is just that in the simple 
Heisenberg model. But for finite Jh, high orders of the 
mediated interaction between core spins make some dif- 
ference. Our result (7) agrees with the conclusions from 
random phase approximation^^, which provides an in- 
tegral equation for the dispersion relation. The strength 
of induced ferromagnetic interaction is determined by 
the hopping energy of the conduction fcrmions. For the 
approach presented to be self-consistent, the spin- wave 
stiffness must be positive. The ferromagnetic order be- 
comes unstable at certain filling when the the stiffness 
vanishes. Though zone boundary spin wave softening 
can be explained by the spin dynamics in the Kondo lat- 
tice model, as shown in this paper. The origin of the 
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